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Abstract 

We consider the chessboard pebbling problem analyzed by Chung, Graham, Morrison 
and Odlyzko [3]. We study the number of reachable configurations G{k) and a related 
double sequence G(k,m). Exact expressions for these are derived, and we then consider 
various asymptotic limits. 
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1 Introduction 

The following problem has attracted some attention recently. We begin with an infinite 
chessboard, which consists of the lattice points '■ i,j > 0} in the first quadrant. We 
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refer to an individual lattice point as a "cell". We start with a single pebble placed at (0,0). 
The first "step" consists of removing this pebble and placing two pebbles at the cells (0, 1) 
and (1, 0). At each subsequent step we remove a pebble at cell and place two pebbles at 
cells (i + 1, j) and + 1), provided that the latter two cells are unoccupied. We consider 
all possible choices of After k steps, there will be a total of k+ 1 pebbles on the board, 

in various arrangements or configurations. 

We let G(k) be the total number of reachable configurations with k pebbles. We define 
the level sets L{1) = : i + j — I}, so that Ui> L(l) is the entire quadrant. The original 

problem, posed by Kontsevich pQ, was to show that L(l) UL(2) UL(3) is unavoidable, in that 
such a set must contain at least one pebble for any reachable configuration. A partial analysis 
of this fact was published thereafter by Khodulev [2]. The first complete proof was given by 
Chung, Graham, Morrison and Odlyzko [3]. It was also shown in [3] that L(l) U L(2) is an 
unavoidable set, as well as certain properties relating to the number of unavoidable sets with 
k pebbles, including the geometric growth rate of this quantity as k — > oo. In jl], Knessl 
obtained some further asymptotic properties relating to the enumeration of unavoidable sets. 
Various extensions of this problem are studied by Eriksson [5] and Warren [6]. 

Suppose we allow more than one pebble per cell and start with an initial configuration of 
one pebble in cells (0, m + 1) and (m + 1, 0) and two pebbles in each of the cells (1, m), (2, m — 
1), (m — 1, 2), (m, 1). Thus there are a total of 2m + 2 pebbles in the level set L(m + 1), 
and we assume that L(M) are empty for M > m + 1. Again the pebble at can only be 
moved if the cells + and (i, j + 1) are empty. Let the number of reachable configurations 
corresponding to this starting arrangement be denoted by G(k, m). In [3], a recurrence relation 
for G(k,m) is derived, and it is shown that, for k > 2, G(k, 0) = G(k). It is also established 
that as k -> oo, G(k) ~ c*a fc where a = 2.321642199494 • ■ ■ and c* = 0.12268707 •■• . These 
constants are characterized in terms of a continued fraction representation of the generating 
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function of G(k). In [TJ, Knessl obtained a more explicit analytic characterization of the 
growth rate a of the number of reachable configurations, and showed that this constant can 
be obtained by solving a transcendental equation that involves two series, each resembling 
Jacobi elliptic functions. Other asymptotic properties of G(k,m) for k and/or m — > oo are 
also established in [7J. 

In this paper, we derive an exact expression for the number of reachable configurations 
G(k) and G(k,m). Using the exact expression, we recover the asymptotic results given in [7J, 
and obtain more explicit analytic expressions for each asymptotic scale. These asymptotic 
scales are k — > oo with m = 0(1); k and m — > oo simultaneously with 2k/m 2 > 1; and 
k, m — > oo with I = k — m(m + 5)/2 = 0(1). 

The paper is organized as follows. In Section 2, we state the basic equations and summarize 
the main results for G(k,m) in Theorem 2.1 and 2.2, and for G(k) in Corollary 2.1 and 2.2. 
In Section 3, we provide brief derivations. 



2 Summary of results 

It is shown in [3] that G(k,m) satisfies the following recurrence equations: 

G{k,0) = 2G(k-l,0) + G(k,l) + 8(k,2) (2.1) 
G(M) = G(k - 3,0) + 2G(k - 2,1) + G(k- 1,2) + G(k- 4,1) (2.2) 
G(k,m) = G(k -m- 2, m- 1) + 2G(k -m- 1, m) + G(k - m, m + 1), m> 2.(2.3) 

Here S(i,j) is the Kronecker delta symbol. As shown in [7J, the boundary condition (12. ip can 
be replaced by 

k 

G{k, 0) = 2 k ~ 2 + 2 fc ^O(/, 1), k > 2. (2.4) 
i=i 
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Using (I2.4p in (12. 2p . G(k,0) can be eliminated, which leads to 

t 

G(k, 1) = 2G(k - 2, 1) + G{k -1,2) + G(k - 4, 1) + 2 k ~ 5 + ^ 2 k - l ~ 3 G(l, 1), k > 5. (2.5) 

i=l 

We introduce a generation function of G(k, m), with 

oo 

V m (z) = {-l) m z- m J2 G ( k > m ) zk - ( 2 - 6 ) 

fc=0 

By explicitly solving for V m (z), we obtain the following integral representation for G(k,m). 

Theorem 2.1 The number of reachable configurations G(k,m) has the following exact ex- 
pression 

G(k, m) = ^- j(-l) m z m - k - l V m {z)dz. (2.7) 

Here C is a closed counterclockwise contour around the origin in the z-plane, with \z\ < 1/2 
on C, 

U+m(m+l)/2 00 m i i 



v 7 n=l L=0 L=l v ' 

and 

S(z) = (2z 2 -3z + 2)5i(z) - (4^ 2 - 4z + 1)5 2 (^) + 2z 2 - z - 1, (2.9) 

where 

oo i 1 

S k (z) = ^ ( _l)m^/2 + (2,-l)V2 "Q _L_^. (2 . 10) 

i=l j=l ^ ^ ' 

Since the total number of reachable configurations for the original problem is G(k) = 
G(k, 0), using (12. 7p in (12.41) with m = 1 we get the exact expression of G{k): 



Corollary 2.1 An exact expression for the total number of reachable configurations G(k) is 

G{k) = 2 k - 2 + — V^dz, 



2iri J c 1 — 2z 

where V\ (z) is given by 112. 8\) with m = 1 . 

Using (12. 7ft . we obtain the following asymptotic expressions for G(k,m). These asymp- 
totic results were first obtained in [7], but there some parameters could only be determined 
numerically 

Theorem 2.2 The number of reachable configurations G(k,m) has the following asymptotic 
expressions: 

1. When k — > oo and m = 0(1), 

m(m+3)/2— k oo m n—1 , 

G(k, m) ~ z * j2(-v n+1 *: {n+1)/2+nm n r~z^ n TT^p' (2 - n) 

6 \ Z *> n=l L=0 1 - Z * L=l y L ~ Z *) 

where z* > is the unique root of S(z) = for \z\ < 1/2 and S'(z*) is the first derivative 
of S(z) evaluated at z = z*. The numerical approximation of z* to 15 decimal places is 
given below 

z* = 0.43072 95931 37930 

2. When k, m — > oo with 2k/m 2 > 1, 

m(m+5)/2— k+1 oo .. 

G(k, m)~ -^—g^Tj II ITTJ+T' P.12) 

Li — U 

3. When l = k- m(m + 5) /2 and 2 < I < m + 3, 

' " 1 °° 1 

Hz) n jz z l+i 



S(z y 

v ' L=0 



(2.13) 
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When k and m — y oo with I > m + 3, A2.13\) holds asymptotically, and then reduces to 

Now we consider the asymptotic expression for G(k) when k — y oo. Since < 1/2, we can 
neglect the term 2 k ~ 2 in the right of (I2.4p and use (12. lip in (12. 4p with m = 1. Since < 1/2, 
we have 

1=1 * 

and hence the following corollary. 

Corollary 2.2 T/ie total number of reachable configurations G(k) has the following asymp- 
totic expression, when k — y oo: 



G{k) 



1 - 2z*)S'(z*) 



oo 



i+l n(n+3)/2 n-_l 



; (l - z;)(i - zl M ) y (l - z^) 2 



i/ere S(z) is as in 112. 9\) and Ii2.10\) . 



3 Brief derivations 

Using the generating function f]2. 6f) in (12. 3p . we obtain the following recurrence equation for 

V m {z) 

V m+1 (z) + V m ^(z) = (2 - z~ m - ^VUz). 

We notice that this equation is of the same form as (3.12) given in [7j, with z = 1/a. Two 
linearly independent solutions are given by (3.29) and (3.30) in [7J. We reject the growing 
solution given by (3.30) in [7j since we expect that V m (z) will be bounded as m — y oo. Hence, 
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we have 



V m {z) 



C(z)—^-exp 
log z 



(m-J) 2 /2-(m-J)/2 



(2 \ 00 
ta 7 J=m+1 
m j 00 1 

n 1 _ n 

L=m+1, 

j ^(-l) n+1 Z -" 2 /2+"/2 JJ £ 

' n=l 



L=0 L=m+1,L^J v 

2 N 00 



L—J\2 



. . 27T / 7T 

<-(^)l eX P 771 

log 2! \ 2 log 2 

1 



■L—n 



X 



n 



L=l,L^n 



(1-Z 



L-n\2 ' 



(3.1) 



where C(z) is independent of m. We note that we can rewrite the two products in (13.11) as 
follows: 



and 



L=0 



n 



L—n 



1 



{-l) m+L z 



m+1 (n+m/2)(m+l) 



L=0 



(l-S 



L—n\2 



Z 



n(n— 1) 



Using (E2D and (E3D in (El]) leads to 



n (i _ z l ? n (1 _ 



ZA2 ' 



(3.2) 



(3.3) 



Kn(^) = C(^)7^-exp / 



log z \2 log 2 / j-J- ( 1 — Z 



00 1 J- 00 

J~J 1 J ^(_l) n + m z n2 /2-"/2+(''+'n/2)(™+l) 



n=l 



X 



m j 
L=0 



rn-1 



n 

L=l v y 



(3.4) 



To determine C(z), we use (I2.6P in the boundary condition ( 12. 5ft . which yields 



1 - 2^ 



1 - 2z 



(3.5) 
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We note that (13. 5p holds for \z\ < 1/2. We introduce the function U m (z) 



n+m n 2 /2-n/2+(n+m/2)(m+l) 

A/ 



n=l 



m | n— 1 | 

II \- Z L+n II (T3>) 
L=0 L=l K ' 



and then V m {z) in H3.4[) is 



V m (z) = C(z)-^—exp 1 



logz \21ogz 



oo 



U m (z) 



(3.6) 



Using (13. 6p with m = 1 and m = 2 in (13. 5p . we obtain C(z) as 



C7(z) = (-2z 5 + z i -3z d + 2z 2 + 2z-l)U 1 (z)-z 2 (l-2z)U 2 (z) 



-i 



log 2 

x exp 

2tt f 



7T 



2 log Z 



\ oo 

yn(i-z^ 

' L=l 



(3.7) 



Instead of using (13. 7p in (13. 6p . we introduce the functions Sk(z) defined in (12.101) and rewrite 
C(z) in terms of the Sk{z). (The reason for making this change is that it allows us to verify 
the equivalence of S(z) = with equation (3.39) in [7], which we will discuss later.) Using 
(I2.10p . U\(z) and U 2 (z) in (13. 7p can be expressed as 



U x (z) = -Sxiz) + (l+ -)s 2 (z) - -S 3 (z) 
V z) z 



(3.8) 



and 



U 2 (z) = -S l (z) + (l + ^ + ^)S 2 (z) 



i + 1 + + Is 4 (s). (3.9) 



We rewrite S k (z) in fl2TT0|) as 



oo i+1 
i+l i 2 /2+(2fe-l)i/2 



s k (z) = j2(-iy-z 

i=l j 
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and then, after some calculation, we obtain the following recurrence equation 



z l - k S k ^{z) + (1 - 2z l - k )S k {z) + z l - k S k+1 {z) = 1. 



(3.10) 



Hence, we can use (13.101) with k = 2 and k = 3 to eliminate £3 (z) and S^z) from (13.81) and 
( 13. 9j) . Thus, after some simplification, we obtain C(z) as 



where S(z) is defined in (12. 9p . Using (I3.1ip in (13.61) . we obtain ( 12 .8p in Theorem 2.1. Then 
we have the integral expression of G(k, m) in ( 12. 7h . 

In the remainder of this section, we discuss the asymptotic approximations to G(k, m). We 
first consider k —} 00 and m = 0(1). By using Rouche's theorem, or by plotting S(z) given 
in ( 12.91) numerically, we notice that there is a single real root for \z\ < 1/2. We denote this 
single root as z*, which satisfies S 1 ^*) = 0. Then in the z-plane with \z\ < 1/2, the integrand 
in ( 12 .7p has a simple pole at z = z*, which is the dominant singularity. We use the residue 
theorem to evaluate the integral in (12.71) asymptotically, which yields (12.111) . We note that 
it's not hard to verify that S(z) = is equivalent to (3.39) in [7j, after we substitute 1/z for a 
in (3.39). Thus z^ 1 = a, whose numerical approximation to 100 decimal places is given in [7]. 
By comparing (12. lip in this paper with formula (3.37) in [7], we obtain an analytic expression 
for the constant C\ in (3.38) in [7] as follows: 



By evaluating (I3.12p numerically, we find that the numerical values of C\ and C\K^ provided 
in (3.42) and (3.43) in [7J are only correct to about 5 decimal places, though they are given 




(3.11) 




(3.12) 
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to 15 places. The correct values to 15 decimal places are 



c x = 2.02740 20474 68498 • • • 

and 

cxK* = 0.28777 77049 35052 

In the limit when k, m — > oo simultaneously, the n = 1 term in ( 12. lip dominates, which 
yields ( 12. 12ft . This recovers the result given in (4.18) in [7], but now with c\ computed 
explicitly. 

Next, we consider k and m large with I = k — m(m + 5)/2 = 0(1). Note that I = 2 
corresponds to the smallest number of possible configurations, as G(k,m) = for I < 1. 
Following [7] we denote G(k,m) as W(l,m) and rewrite ( 12.71) as 



n (n-l)(n/2+m+l) 



m | n—1 



For a sufficiently small closed contour C around the origin in the z-plane and I > 2, z = is 
the only pole inside of C, and it is of order I — 1. Thus, using the residue theorem in (13.131) 
leads to 



^ = (r^£5^{s 5 )E(-i)^"- 1,( " /2+ " ,+1) 



(Z-2)! *->o 1 5^) ^ 

/ v \ / n= x 

n n 713^ • 
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The first two terms from the infinite sum in the right-hand side of (I3.14p are 



J z m+2 m ^ 

n i _ z l+i + a _ n 



L=0 



L+2 ' 



(3.15) 



L=0 



and the remaining terms are of order 0(z 2m+5 ) as z — > 0. We notice that after taking / — 2 
derivatives, the second term in (13. 15[) is of order 0(z m ~ l+A ) as z — > 0. This implies that as 
long as m — I + 4 > 1, i.e., I < m + 3, only the n = 1 term in the infinite sum in f)3.14p 
contributes to the derivative at z — 0, which yields 



lim 



1-2 



(Z-2)!^o cZz'- 2 [5(2) 



1 m 1 

Yi) n i - z l+i 



2 < I < m + 3. 



(3.16) 



Rewriting the expression in the brackets in (I3.16P as 



I in | 11 x 

sTz) n i _ z l+i = s(zj n i _ z l+i n ^ - 

v y L=0 v 1 L=0 L=m+l 



we see that (13.161) is independent of m. This leads to (I2.13p . and gives an analytic expression 
for Wo (I), which appears in [7] but there no explicit expression is given. This concludes our 
derivation. 
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